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We study the effects of insulating oxides in their crystalline forms on the energy band structure 
of monolayer and bilayer graphene using a first principles density functional theory based electronic 
structure method and a local density approximation. We consider the dielectric oxides, Si02 («- 
quartz) and AI2O3 (alumina or a-sapphire) each with two surface terminations. Our study suggests 
that atomic relaxations and resulting equilibrium separations play a critical role in perturbing the 
linear band structure of graphene in contrast to the less critical role played by dangling bonds that 
result from cleaving the crystal in a particular direction. We also see that with the addition of a 
second graphene layer, the Dirac cone is restored for the quartz surface terminations. Alumina needs 
more than two graphene layers to preserve the Dirac cone. Our results are at best semi- quantitative 
for the common amorphous forms of the oxides considered. However, crystalline oxides for which 
our results are quantitative provide an interesting option for graphene based electronics, particu- 
larly in light of recent experiments on graphene with crystalline dielectrics (hexagonal BN) that 
find considerable improvements in transport properties as compared to the those with amorphous 
dielectrics. 

PACS numbers: 71.15.Mb, 71.20-b, 73.20.-r 



I. INTRODUCTION 



Recent advances in large area graphene films on 
metal substrates and understanding the morphology of 
its domains^, deposition of dielectric^, making metal 
contacts^ and designing novel electronic switches^ are 
making graphene based electronics more likely. Theo- 
retical efforts have played a crucial role in this devel- 
opment often elucidating the complex interplay of sub- 
strates and dielectrics with graphene. Notable among 
these efforts is the understanding of the role of interface 
and substrate charge inhomogeneity in degrading the car- 
rier mobilities^. Density functional theory (DFT) based 
electronic structure methods, used in some of these stud- 
ies, have provided inputs for optimization of various pro- 
cesses used in graphene electronics. A crucial obstacle 
in graphene electronics research has been the need to im- 
prove the substrate quality to maintain the novel intrinsic 
graphene properties. In this article, we address the inter- 
face electronic structure of graphene with two crystalline 
oxides namely, quartz (Si02) and alumina (AI2O3), us- 
ing ab-initio DFT based numerical method. Crystalline 
substrates/dielectrics provide a compelling alternative to 
amorphous substrates/dielectrics as they can provide an 
order of magnitude improvement in the electronic trans- 
port properties compared to the the amorphous oxides^. 
Our studies have implications in interpreting future ex- 
periments with crystalline oxides. 

The atomic structure of graphene on amorphous in- 
sulating substrates has been explored in some experi- 
mental studies^ using nanometer scale microscopic tech- 
niques such as scanning tunneling microscopy and atomic 
force microscopy. These studies suggest that graphene 
conforms to the structural inhomogeneity of the un- 
derlying substrate. In addition, a study using Raman 
spectroscopy examined the role of amorphous oxides on 



phonon dispersions, indirectly hinting at the change in 
the electronic spectrum of graphene due to the presence 
of oxide substrates^. We are aware of two recent DFT 
studies of monolayer graphene on a crystalline Si02^»5^, 
which discuss perturbations to the linear band structure 
due to the presence of oxides, the role of dangling bonds 
and their passivation in influencing the linear band struc- 
ture. 




FIG. 1. (Color online) Schematic illustrations of supercell 
structure of monolayer graphene on crystalline oxides, (a) 
Graphene on Si-terminated quartz with four unit-cells, (b) 
Graphene on Al-terminated alumina with four unit-cells. The 
atoms are shown in color: Si(blue), Al(purple), O(red) and 
C (yellow). 



Our paper is organized as follows. In section II, we dis- 
cuss the computational method and convergence param- 
eters used for this study, followed by the motivation and 
procedure used to build surface models from their bulk 
counterparts. The results from these calculations are dis- 
cussed in section III, which details the role of atomic 
relaxations, the resulting equilibrium distances and the 
role of unsaturated or dangling bonds of oxide surfaces in 
perturbing the linear band structure of graphene. We ad- 
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dress the origin of change in hnear band dispersions with 
the help of orbital and atom projected densities of states 
and also present the energetics of our interface models. 
Finally, we summarize our results and present our con- 
clusions. 



II. COMPUTATIONAL METHOD AND 
SURFACE MODELS 

This section addresses the details of the compu- 
tational method we used, followed by the procedure 
adopted to obtain the surface models of quartz and 
alumina from the bulk and the convergence parame- 
ters employed. We adopted a plane-wave based elec- 
tronic structure method^ ^ using a local density approx- 
imation (LDA)^^ for exchange and correlation, and the 
projector augmented plane-wave potential for electron- 
ion interaction^^. The bulk structures for both quartz 
and alumina are consistent with those mentioned in 
literature^^ and the surface models built from the bulk 
structures conform to the widely accepted quartz and 
alumina structured containing alternating cation and 
anion layers. 

The unit cell of bulk quartz contains 27 atoms with 3 
silicon planes and 6 oxygen planes, each plane comprising 
three atoms. The unit cell of alumina contains 30 atoms 
with 4 aluminium planes and 6 oxygen planes, again with 
each plane consisting of only three atoms. We used a 7 
X 7 X 5 k-point mesh in the hexagonal Brillouin zone 
(BZ) and a kinetic energy cut-off of 612 eV. The results 
were carefully checked in comparision to a larger k-point 
set and higher energy cut-offs. 

After optimization of atomic positions and lattice con- 
stant for the bulk structures, we found these values to 
be very close to the experimental values, differing by less 
than 1% (Table I). Using these optimized lattice param- 
eters, we constructed surface models of both quartz and 
alumina. Four oxide bulk unit cells were chosen to rep- 
resent the thin-films and were stacked along the crystal- 
lographic c-direction (corresponding to a film thickness 
of about 22 A). At their top termination, we placed a 
graphene layer. We found that 6 x dc-c graphene, con- 
taining 24 carbon atoms (where dc-c = 1-42 A), was 
nearly commensurate with the hexagonal surface of the 
oxides (Fig. 1). With this set-up, we get the lattice mis- 
match of quartz and alumina with graphene as ^ 0.19% 
and ~ 0.42%, respectively. The lattice mismatch values 
mentioned here serve only as an initial guideline towards 
assessing the degree of distortion of graphene in presence 
of the oxides. We keep these mismatched values of lat- 
tice constants by not allowing any variation in supercell 
size. These mismatched values result in a weak rippling 
effect of carbon atoms which may not be observed in ex- 
periments using graphene on atomically flat substrates. 
We estimate average height fluctuations of relaxed two- 
dimensional graphene and find that the deviations are 
approximately 0.05 A, which is small compared to fluc- 
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FIG. 2. Energy band structures of (a) monolayer graphene 
on Si-terminated quartz (b) monolayer graphene on O- 
terminated quartz (c) bilayer graphene on 0-terminated 
quartz. The atoms in the supercell were relaxed. The Fermi 
energy is set at zero. The origin of occurrence of linear bands 
at the T-point instead of i^^-point of the supercell is explained 
in the text. 



tuations in suspended graphene^. 

Modehng roughness present on amorphous oxide sur- 
faces by DFT is challenging partly due to the inherent 
complexity in generating the amorphous structures by 
DFT (the traditional melt and quench method) as well as 
the requirement of relatively large size of the supercell to 
represent roughness, a computationally intensive task for 
DFT calculations. Moreover, the amorphous oxide struc- 
tures do not possess dangling bonds. Therfore, an addi- 
tional constraint in using the crystalline substrates, be- 
sides treating surface roughness in an approximate way, 
is to saturate these dangling states for assessing the ex- 
ternal effects on graphene electronic structure. To make 
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the combined graphene and oxide thin films structure 
an isolated system, we used a vacuum region adjacent 
to graphene layer in the stacking direction and supercell 
is then repeated in all three crystallographic direction 
by imposing periodic boundary conditions on the crys- 
tal wave functions. The dangling states at the bottom 
of the supercell were saturated with hydrogen atoms and 
the atoms within the top two unit-cells of the oxide film 
and all the graphene atoms were allowed to relax while 
keeping the cell size fixed. We used the same energy cut- 
off value as in the bulk calculations, but the k-point mesh 
in the BZ was chosen to be 7 x 7 x 1. The convergence of 
results with respect to larger k-mesh size, larger vacuum 
sizes and energy cut-off was also checked. For atomic re- 
laxation, the total energy was assumed to have converged 
when all the components of the Hellman-Feynman forces 
were smaller than 0.01 eV/A. 




Wave-vector 

FIG. 3. Energy band structures of (a) monolayer graphene 
on Al-terminated alumina (b) monolayer graphene on O- 
terminated alumina. The atoms in the supercell were relaxed. 
The Fermi energy is set at zero. 



III. ENERGY BAND DISPERSIONS, 
DENSITIES OF STATES AND ENERGETICS 

In this section, we discuss the effect of underlying ox- 
ide surfaces on the electronic band structure of graphene 
for different surface terminations. In particular we ex- 
plore the impact of graphene-oxide distances, oxide dan- 
ghng bonds and atomic relaxations on the band struc- 



ture. With relaxations of atomic positions, we get equi- 
librium oxide-graphene distances tabulated in Table I. 
All the distances are larger than the corresponding in- 
teratomic distances of the constitutent atoms in bulk 
or molecular compound phases, suggesting that forma- 
tion of compound phases at low temperature is unlikely. 
However, higher temperture effects may support forma- 
tion of these phases. We present detailed examination 
for all four surface terminations, beginning with silicon- 
terminated quartz. Figure 2(a) shows the self-consistent 
band structure of graphene on silicon-terminated quartz 
at the equilibrium distance (Table I). It is evident that 
hnear band structure is not perturbed at such a distance, 
which is consistent with atom and orbital projected den- 
sity of states (DOS) (figure not shown). We find weak hy- 
bridization between carbon-p and silicon-p orbitals in the 
valence band region up to the Fermi level. We note here 
that the Dirac point occurs at the center of the supercell 
BZ and not at the K point of the honeycomb BZ due 
to symmetry reasons explained in the Appendi:?^^^. For 
silicon-terminated quartz, the graphene sheet remains 
flat. On oxygen-terminated quartz, our calculations sug- 
gest noticeable rippling of graphene. Moreover, the small 
equilibrium interface distance of 1.76 A (Table I) leads 
to perturbation in graphene's band structure (Fig. 2(b)). 
The DOS plot (Fig. 4(a)) shows considerable hybridiza- 
tion in the whole valence band region between oxygen-p 
and carbon-p orbitals. A similar conclusion was reached 
in a recent DFT-based calculation of monolayer graphene 
on oxygen-terminated quartz Interestingly, adding one 
more layer of graphene to the oxygen-graphene supercell 
allows the first layer to act as a buffer, protecting the 
Dirac cone (Fig. 2(c)). Both surface terminations of alu- 
mina result in strong perturbations of the graphene Dirac 
cone (Fig. 3(a) and (b)) which is consistent with the DOS 
plots (Fig. 4(b) and 4(c)). We note that adding an extra 
layer of graphene to the alumina surface terminations 
does not recover the linear band structure, suggesting 
more than two graphene layers are needed to protect the 
Dirac cone. 

To understand the distance-dependent electronic 
structure of the combined graphene and oxide system, 
we manually adjusted the interface distances both be- 
low and above the calculated equilibrium distances. The 
atomic positions were not relaxed while other compu- 
tational parameters kept the same as those used in the 
self-consistent runs discussed above. We demonstrate our 
calculations by using silicon-terminated quartz as an ex- 
ample. We chose two representative distances 2 A and 4 
A which are, respectively, below and above the equilib- 
rium distance of 3 A (Table I). The Dirac cone is pro- 
tected in 4 A case (Fig. 5(a)) but not in 2 A case (Fig. 
5(b)), suggesting a strong dependence of band structure 
on the interfacial distance. We have performed these 
tests with other surface terminations and found that the 
interfacial distance dictates whether the Dirac cone can 
be protected or not. The equilibrium distance after relax- 
ation is thus quite important in determining graphene's 
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properties. 
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FIG. 4. (Color online) Atom and orbital projected den- 
sity of states for a C atom in monolayer graphene and the 
atoms of the topmost planes in (a) 0-terminated quartz (b) 
Al-terminated alumina and (c) 0-terminated alumina oxides. 
The atoms in the supercell were relaxed. The Fermi energy 
is set at zero. 

Table I shows the binding energy values for graphene 
on four surface terminations considered in this study. In 
the most relevant case, i.e., no surface passivation but 
with atomic relaxation, these values hint at the non- 
bonding nature of graphene to the underlying oxide ox- 
ides. The binding energy values are obtained by using 
the following definition. 



E 



bind 



E(supercell) - E{Gt) - E(oxide) (1) 



where E'(supercell) denotes the total energy of the su- 
percell containing the oxide and a graphene layer. E'(Gr) 
and oxide) denote, respectively, the total energies of 



isolated graphene and isolated oxide in the same supercell 
set-up, with the same energy cut-off and k-point mesh as 
that of the combined graphene and oxide calculations. 



TABLE I. Lattice parameters (in A) of bulk quartz and alu- 
mina, average interplanar distances (in A) of graphene from 
the four underlying surface terminations and their binding en- 
ergies (in eV/atom), in case of non-passivated but relaxed sur- 
faces. The numbers in parenthesis are the out-of-plane lattice 
constants of the bulk phases, and those in square parenthesis 
are interatomic distances between carbon and the correspond- 
ing atom of the surface terminations in the bulk or molecular 
phase. It should be noted that for isolated oxide calculations 
we passivated the top layer with hydrogen atoms to keep the 
supercell non-magnetic. 



Lattice Paramters 



d(C-x) 
(x = Si,Al,0) 



Bulk quartz 
This work 
Expt.^ 
Si-terminated 
0-terminated 



4.914 (5.408) 
4.913 (5.405) 



3.0 (1.89)^ 
1.76(1.3)5 



10.061 
0.581 



Bulk alumina 
This work 
Expt.^ 
Al-terminated 
0-terminated 



4.907 (4.908) 
4.943 (4.907) 



2.7 (1.89-2.19)^ 
2.15 



1.017 
0.689 



^ Reference 17, 
^ Reference 22 
^ Reference 23 
^ Reference 18, 
^ Reference 25 



21 



24 



To assess the importance of danghng bond saturations 
on graphene's electronic specturm, we performed a calcu- 
lation with the top oxide surface saturated with hydrogen 
atoms. This calculation was undertaken only for oxygen- 
terminated surface of quartz since silicon-terminated sur- 
face already protects the Dirac cone without surface pas- 
sivation. Both aluminium- and oxygen-terminated sur- 
faces of alumina, even in the initial structures, show non- 
trivial stoichiometry, as a result passivation of dangling 
bonds is avoided. Two separate calculations were per- 
formed with equilibrium oxygen-carbon interfacial dis- 
tance in oxygen-terminated quartz: one with un-relaxed 
and unpassivated oxide surfaces, and another with un- 
relaxed passivated oxide surface with graphene on top of 
it. We find that on un-passivated surfaces, there is no ef- 
fect on graphene's band structure except for some doping 
effects (Fig. 6). This hints at the non-critical role played 
by dangling bond states. To check the effect of relax- 
ation of the unpassivated surfaces in all four terminations 
without the graphene layer, we performed additional cal- 
culations which show that both aluminum- and silicon- 
terminated surfaces are inherently stable compared to 
the initial surface model we built from the correspond- 
ing bulk structures as well as the final relaxed unpassi- 
vated structures. The interatomic and interplanar dis- 
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FIG. 5. Energy band structures of (a) monolayer graphene 
on Si- terminated quartz at a distace of 4 A (b) monolayer 
graphene on Si-terminated quartz at a distance of 2 A. The 
atoms in the supercell were not relaxed. The Fermi energy is 
set at zero. The oxide surface was unpassivated. 
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FIG. 6. Energy band structures of (a) monolayer graphene on 
0-terminated quartz at a distace of 1.76 A. The top surface 
of oxide is passivated with hydrogen atoms (b) monolayer 
graphene on O- terminated quartz at a distance of 1.76 A. 
The top surface of oxide is not passivated in this case. The 
atoms in the supercell were not relaxed. The Fermi energy is 
set at zero. 



IV. SUMMARY AND CONCLUSION 



tances were found to be not different from the unre- 
laxed initial and relaxed final structures and the surface 
atomic configurations were similar. This situation corre- 
sponds to the configuration adopted by the silicon and 
aluminium surface terminations in experiments immedi- 
ately after cleaving from the bulk crystals. Therefore, 
our initial surface models with graphene on the top of 
these surfaces and the corresponding interfacial distances 
in Table I obtained upon relaxing these combined struc- 
tures without performing the intermediate step relaxing 
the unpassivated oxide surfaces, are valid. The situation 
is not same with oxygen-terminated surfaces of quartz 
and alumina. Both surfaces show large relaxations of 
in-plane surface atoms. However, they are all displaced 
by the same amount. To a first order approximation, 
this uniform displacement is equivalent to displacing the 
graphene layer in-plane with respect to the oxide sur- 
face. Since our graphene position is an independent vari- 
able, this displacement would not matter. We believe 
all second-order effects of relative in-plane movement be- 
tween substrate layers are taken care of in the relaxations 
of the combined structures that we perform. 



In summary, we have studied, using a first principles 
DFT method, the effect of two crystalline insulating ox- 
ides, quartz and alumina, on the electronic structure of 
monolayer graphene. Silicon-terminated quartz retains 
the graphene band structure on non-passivated surface 
with atomic relaxations. Oxygen- terminated quartz and 
both alumina terminations perturb the linear bands of 
graphene. Graphene and oxide interface distance is criti- 
cal in determining the extent of disturbance to graphene. 
By comparison, the effect of dangling bond states on 
the oxide surface is not very crucial. Bilayer graphene 
on oxygen-terminated quartz retrieves its linear band 
structure while graphene on alumina terminations (Al 
or O) perhaps needs additional graphene layers. Our 
studies are relevant in situations where a crystalline di- 
electric is used either as a template for growth or as a 
substrate/dielectric for exfoliated graphene experiments. 
Recent use of crystalline BN as a substrate/dielectric for 
graphene experiments is a promising step in this direc- 
tion, and we believe more experiments with crystalline 
oxides should be conducted in future, which make DFT- 
based studies useful for interpreting experiments. 
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Appendix 

Quartz, alumina and the graphene layer, all have a 
hexagonal unit cell. Let Ai, A2 and As be the lattice 
vectors of the entire graphene and oxide system and 5i, 
^2, and Bs be the correponding reciprocal lattice vec- 
tors. Similarily, let ("ai, "a 2 , ^3) and (61, 62 and 63) 
be the triad of primitive and reciprocal lattice vectors of 
the graphene layer alone. 

We note that the lattice structure of the graphene 
plus oxide system is a (2 \/(3) x 2 y^3))R30^ recon- 
struction of the lattice structure of the graphene layer 
alone. Thus A = 2y^3)a where A is the lattice con- 
stant of the entire system and a is the lattice constant 
of the two-dimensional graphene layer. Since both unit 
cells are hexagonal this implies that the reciprocal lat- 
tice structure of the graphene layer on top is a (2y^3) x 
2 Y^3))R30^ reconstruction of the reciprocal lattice struc- 
ture of the graphene plus oxide system below. Taking "ai 
to be along x. 



ai = ax 



^2 = ail±+^y) 



Ai = 2^/(3)a(4^x+iy) 



A2 = 2\/(3) 



As = cz 



ay 



(A.1) 
(A.2) 

(A.3) 

(A.4) 
(A.5) 



oi = — (x F — y) 



^2 = — (^=::7y) 



« V(3)^ 




(A.6) 



(A.7) 



(A.8) 



FIG. 7. (Color online) Schematics of supercell Brillouin zone 
(small) superposed on the graphene honeycomb Brillouin zone 
(large). The reciprocal lattice vectors are shown along with 
the overlap of F to K vector of graphene with the Bi vector 
of the supercell. The overlap is due to the symmetry reasons 
and is explained in the text. 



B. = 



27r 



2\/(3)a 7(3) 



-x + y) 



(A.9) 



The reciprocal lattice vectors are. 



Figure 7 shows the BZ set-up with these vectors using 
equations A.6 to A.9. 

Now if we take the same F point as the center of both 
the Brillouin zones and draw the reciprocal unit cell of 
both the structures, we realize that the K point of only 
the graphene layer lies on the reciprocal lattice vector Bi 
of the graphene plus oxide system (as shown in figure). 
This means that the K point of graphene layer folds in 
by symmetry onto the F point of the entire system. 
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